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Over the last years both cosmic-ray antiproton measurements and direct dark matter searches 
have proved particularly effective in constraining the nature of dark matter candidates. The present 
work focusses on these two types of constraints in a minimal framework which features a Majorana 
fermion as the dark matter particle and a scalar that mediates the coupling to quarks. Considering 
a wide range of coupling schemes, we derive antiproton and direct detection constraints using the 
latest data and paying close attention to astrophysical and nuclear uncertainties. Both signals are 
strongly enhanced in the presence of degenerate dark matter and scalar masses, but we show that 
the effect is especially dramatic in direct detection. Accordingly, the latest direct detection limits 
take the lead over antiprotons. We find that antiproton and direct detection data set stringent lower 
limits on the mass splitting, reaching 19% at a 300 GeV dark matter mass for a unity coupling. 
Interestingly, these limits are orthogonal to ongoing collider searches at the Large Hadron Collider, 
making it feasible to close in on degenerate dark matter scenarios within the next years. 



I. INTRODUCTION 



Complementarity between different signals has always 
been highlighted as the preferred way forward in the 
search for weakly interacting massive particles (WIMPs) 
- see [IHI] for general-purpose reviews. However, WIMP 
searches were for a long time an extremely data-starved 
field of research with experiments lagging far behind the- 
oretical predictions. That is no longer the case. With 
large amounts of data pouring in from direct, indirect 
and collider dark matter (DM) searches, complementar- 
ity studies are now possible at an unprecedented level of 
detail [MS]. Direct searches, in particular, have devel- 
oped tremendously over the last decade benefiting from 
the use of various targets and techniques to measure 
WIMP-induced nuclear recoils. Presently, the situation is 
rather blurry: while experiments such as DAMA/LIBRA 
[TDHE], CoGeNT [H H and CRESST [Tg report ex- 
cess of nuclear recoil events and also evidence for annual 
modulation in the case of DAMA/LIBRA and CoGeNT, 
other collaborations - including XenonlO/100 [ISlU9| . 
CDMS [Mil], SIMPLE [23HH] or COUPP [251 rj7] - 
find null results. Leaving aside this controversy, that has 
been discussed at length in the literature [28 33 , it is 
clear that the whole array of different targets used boasts 
huge sensitivities to both spin-dependent (SD) and spin- 
independent (SI) WIMP-nucleus scattering. Together 
with direct detection, indirect searches via antiprotons 
have provided valuable hints in shaping the phenomenol- 
ogy of viable WIMP models. In fact, the exquisite data 
on the cosmic-ray antiproton flux collected by experi- 
ments such as BESS [SH36] or PAMELA [37] fall nicely 
on top of the expectations from cosmic-ray spallations 
in the Galaxy. This means in practice that dark matter 
annihilations or decays cannot yield copious fluxes of an- 
tiprotons, and correspondingly the coupling to quarks is 
much constrained [33j E2] ■ Notice that this is precisely 



the coupling that drives WIMP-nucleus scattering in un- 
derground detectors, which makes antiprotons and direct 
searches particularly suitable for complementarity stud- 
ies. 

Several complications arise when exploring the com- 
plementarity between two or more observables. Firstly, 
if we are to draw sound conclusions, all relevant uncer- 
tainties must be treated carefully. In the case at hand, 
the unknowns are sizable and different in nature: on the 
one hand, antiproton searches are prone to uncertain- 
ties on the galactic dark matter profile j40T[43] and on 
cosmic-ray propagation [4"1H5U] ; on the other hand, di- 
rect searches suffer from the lack of knowledge on the 
local dark matter density [5TH5U] and local velocity dis- 
tribution [61-67 as well as from nuclear uncertainties 
[6"3l IBT)] . Secondly, combining signals such as antiprotons 
and direct detection requires the specification of an un- 
derlying particle physics framework. For instance, the 
author of references [70l [71] points out the usefulness of 
antiproton measurements to constrain light WIMPs able 
to accommodate the hints of signal in DAMA/LIBRA 
and CoGeNT. This claim is, though, dependent on how 
dark matter couples to quarks, especially the light quarks 
abundant in nucleons that get struck in direct detection 
experiments. A definitive conclusion is therefore neces- 
sarily model-dependent. This highlights the difficulty in 
pursuing complementarity studies and deriving conclu- 
sive results. 

Now, model building imposes hardly any boundaries 
on the phenomenology of dark matter candidates. Nev- 
ertheless, one can learn a lot by focussing on relatively 
simple realisations. For example, it has been recently em- 
phasized [72H75] that the presence of mediator particles 
degenerate in mass with the WIMP leads to enhanced 
signals in both antiprotons and direct detection. In the 
case of the former, degeneracy boosts the contribution 
of 2 — > 3 processes to the annihilation yields [73H73] . 
while in the latter it induces almost-resonant scattering 
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in SD and SI interactions [75]. Besides, all this proceeds 
at mass degeneracies below the trigger of searches at the 
Large Hadron Collider (LHC). To the best of our knowl- 
edge, the complementarity between antiprotons and di- 
rect searches has not been conveniently explored in this 
important case of degenerate mass states. The present 
work is precisely devoted to fill that gap, and is a first 
step in exploring realistically the interplay between di- 
rect and indirect searches within mass-degenerate WIMP 
frameworks. In particular, we show explicitly how data 
already at hand effectively constrain the presence of de- 
generate mass states, and how a reasonable experimental 
effort over the coming years will be able to close in on 
these scenarios. As we shall see, collider searches are or- 
thogonal to antiproton and direct detection, and hence 
a multi-disciplinary approach is highly desirable for the 
near future of dark matter searches. 

We start by introducing the main ingredients of our 
particle physics framework in Section [II] Then, Sections 
|III| and |IV| are devoted to the formalism behind antipro- 
ton and direct detection constraints, respectively. Tak- 
ing special care in modelling all relevant uncertainties, 
we derive in Section [V] the limits imposed by the latest 
SD and SI direct detection searches as well as cosmic-ray 
antiproton data, before concluding in Section |VI| 



II. PARTICLE PHYSICS MODEL 

Direct detection experiments probe mainly the inter- 
action of dark matter with the light quarks, the funda- 
mental scattering process being xi ~^ XI- As is well- 
known, the corresponding annihilation process XX QQ 
that is obtained from crossing the scattering diagram 
is strongly suppressed for Majorana fermions, either by 
the quark mass or the velocity v ~ 10~ 3 c in the Galac- 
tic halo. Therefore, at first sight, it may seem hope- 
less to constrain the coupling of Majorana dark mat- 
ter to light quarks via indirect detection. However, the 
helicity suppression can be lifted by the additional ra- 
diation of a spin-one boson in the final state [TSj [77] . 
Because of this property, the two-to-three annihilation 
channels XX iffl ■, with V being a gluon [7THT5] , a 
Z/W-boson [7J El IBDHUI (see also [53 [SB]) or a pho- 
ton ESI EH IH7J typically dominate over the two-to-two 
process, provided that the particle mediating the annihi- 
lation has a mass that is of the same order as the dark 
matter mass |73j . In particular, the emission of gluons 
leads to a rather efficient production of antiprotons, mak- 
ing it possible to obtain limits on dark matter coupling 
to light quarks from indirect detection, and consequently 
to relate direct and indirect detection constraints in an 
unambiguous manner. 

In the following, we consider a simple model for mass- 
degenerate dark matter that is suitable for a comparative 
analysis of direct and indirect detection constraints. It 
consists of a minimal extension of the Standard Model 
(SM) by a Majorana fermion x, which constitutes the 



dark matter in the Universe, and a scalar particle r\ 
not much heavier than x, which provides the portal to 
the Standard Model via a Yukawa interaction with the 
light quarks El- Although this simplified model is not 
uniquely related to a particular framework of physics be- 
yond the SM, its particle content is motivated by the min- 
imal supersymmetric extension of the Standard Model 
(MSSM). 

A. Outline of the model 

The dark matter particle x IS taken to be a singlet un- 
der the Standard Model gauge group, while r) is a singlet 
under SU(2), a triplet under 577(3) and carries hyper- 
charge such that the electric charge of the quarks is ob- 
tained. The total Lagrangian of the model can be written 
as: 

£ = C SM + \x c i$X ~ \rn x x c X + {D^(D^) 

- mffir] + L int . (1) 

Here, Lsm denotes the Standard Model Lagrangian, and 
Dp is the usual covariant derivative. The remaining part 
of the Lagrangian, Li n t, contains the interactions of the 
new particles with the quarks. In the following, we as- 
sume that r\ carries a flavour quantum number, thus en- 
suring that r\ couples only to one quark flavour and does 
not induce any flavour changing effects. To be more pre- 
cise we consider the cases where rj has U = — 1, D = — 1, 
S = — 1 or B = — I as well as a model with three mass 
degenerate scalars 7? llj d,s, each of which couples to one of 
the light quarks. Thus the interaction Lagrangian can be 
expressed as 

L m t = -fx^RV + h-c. , (2) 

where VP = u, d, s, or b. Similarly, one could also consider 
an analogous model where dark matter couples to the 
left-handed quarks. However, we do not discuss this case 
separately here since it yields very similar constraints. 

Formally, the simplified model described above can 
be obtained within the MSSM in the limiting case of 
a bino as the lightest supersymmetric particle, and a 
right-handed squark as the next-to-lightest one. The 
Yukawa coupling / = fsuSY is then fixed in terms of the 
U{1) Y gauge coupling g' to fsusr = 4g'/(3\/2) ~ 0.33 
for up-type quarks and Jsusy = 2<//(3v2) ~ 0.16 for 
down-type quarks (for the analogous model involving 
left-handed quarks, one would have fsuSY = sV(3v2))- 
In order to explore the full available parameter space, 
we shall consider the coupling / as a free parameter in 
the following. The quantum numbers of the scalar par- 
ticle allow an additional renormalizable coupling to the 
Higgs field, L\ nt — X^W Hrfr/. It provides an additional 
mass term for r/, which can be absorbed into the tree- 
level mass, mi + ^3 v ew ~~ ^ m n- Apart from that, it 
gives rise to couplings of rj to the Higgs boson, however 
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since only the couplings of the dark matter particle x are 
relevant for direct and indirect detection this turns out 
to be irrelevant for these limits. The most fundamen- 
tal constraint on light coloured particles comes from the 
non-observation of an excess in the invisible decay width 
of the Z boson at LEP, Ar in „ < 2.0 MeV [88 . This 
excludes exotic charged particles with a mass below 40 
GeV [HH]j and since we are interested in the case of a 
compressed mass spectrum we only consider dark matter 
masses m x > 40 GeV throughout this analysis. 



B. Thermal freeze-out 

One of the main arguments in favor of a WIMP as a 
dark matter candidate is that such a particle can be pro- 
duced quite naturally by thermal freeze-out in the early 
Universe. The model and the parameter space of interest 
here, i.e. a compressed mass spectrum of the dark mat- 
ter and the next-to-lightest particle beyond the Standard 
Model, make it necessary to treat thermal production 
with special care, since coannihilation can induce great 
deviations of the relic abundance from generic expecta- 
tions [30] • In order to include all the relevant interactions 
of r\ and \ we choose to use a fully numerical calculation 
of thermal freeze-out by MicrOMEGAS [HT] which we 
have checked against a semi-analytic approximation |92j . 
The analytic result for the freeze-out density can be ex- 
pressed as 
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where N c = 3 is a colour factor, and gqv as well as 
TfQ ~ m x /25 correspond to the annihilation cross-section 
and freeze-out temperature without taking coannihila- 
tions into account [731 [71], respectively. The effective 
cross-section is given by 
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while the cross-sections for the processes XV ~ * 19 y VV ~~ * 
qq, and r]fj — > gg have a non-zero s-wave contribution for 

m„ — > 0, 
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The s-wave contribution for r\f\ — > qq (see Fig. [T] diagram 
(h) as well as (d) with flipped charge flow in one of the 
lines) is helicity-suppressed, and therefore less important. 
For any given mass spectrum, we fix the coupling / such 
that the relic density matches the WMAP value [53] . For 
a large mass splitting, m n — m x 3> TfQ, coannihilations 
are suppressed, and the resulting coupling / is of order 
1/3 for dark matter masses in the 100 GeV range. For 
smaller mass splittings, additional coannihilation chan- 
nels contribute, and therefore a smaller coupling / is suf- 
ficient to achieve the same relic density. Below a certain 
value for the dark matter mass and the mass splitting, 
the contribution of the annihilation channel rjfj — > gg, 
whose cross-section depends only on the strong coupling, 
can become so large that the freeze-out density falls be- 
low the WMAP value even for / = 0. For a mass split- 
ting m v /m x = 1.1 (1.01), this leads to a lower limit on 
the dark matter mass of m x > O[200 (1000) GeV], for 
which thermal freeze-out can account for the observed 
dark matter density. This also explains the behaviour 
of the thermal relic density constraint shown in Fig. [2] 
below. The results of the numerical solution and the ap- 
proximation are in good agreement, thus indicating that 
both approaches yield a reasonable description of ther- 
mal freeze-out. We have also checked that the inclusion 
of diagrams similar to the ones shown in Fig. [T] but with 
gluons replaced by electroweak gauge bosons or Higgs 
bosons, leads to minor corrections, assuming | A3 1 < 1 for 
the coupling to the Higgs. We also checked that anni- 
hilation into top quarks and into a pair of Higgs bosons 
has no sizeable impact. In the following we prefer to 
use the limits obtained by MicrOMEGAS since they are 
expected to be more precise. 



with n1 q = g l (m i T/(2Tr)) 3 / 2 e- m ^ T and g x = 2, g n = 
g r - t = 1. When expanding av = a + bv 2 into s- and 
p-wave contributions, the thermal average is given by 
{av) = a + 3bT/m x . The Feynman diagrams for the 
leading processes which we take into account are shown 
in Fig. [l] The cross-section for the annihilation of dark 
matter particles, XX ^ is helicity- and/or velocity- 
suppressed, 
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III. ANTIPROTON CONSTRAINTS 

In this section we first discuss the annihilation pro- 
cesses that are relevant for indirect detection, and then 
review constraints arising from the primary flux of an- 
tiprotons induced by the decay of annihilation products. 
As discussed above, the lowest order annihilation pro- 
cess XX ™>> QQ? see Eq. ([5]), is strongly suppressed. In 
particular, since the typical dark matter velocity in the 
halo of the Milky Way is of the order v ~ 10~ 3 c, the p- 
wave contribution to the cross-section is suppressed much 
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FIG. 1: Feynman diagrams showing the leading (co) annihilation processes that contribute to the thermal freeze-out density. 
Two additional diagrams corresponding to (a) with reversed charge flow and to (d) with flipped legs are not shown, as well as 
several diagrams for charge conjugated processes. 



more strongly than during thermal freeze-out. On the 
other hand, the s-wave contribution is suppressed by the 
masses of the light quarks, and therefore even smaller. 
Hence, 2 — > 3 annihilation processes, where the quark 
pair is accompanied by a gauge boson, typically yield 
the dominating contribution to the annihilation within 
our Galaxy. The differential cross-section for these pro- 
cesses is given by [71] 



A. Formalism &z uncertainties 

The annihilation of dark matter particles in the halo 
of the Milky Way induces a primary flux of antiprotons 
that is potentially observable at the Earth. The rate of 
antiprotons produced at position r with respect to the 
center of our Galaxy per unit of kinetic energy and vol- 
ume is given by 



vdajxx -> qqV) = Cyf* 
dE v dE q 8tt 2 to4 

(Xmax -x)[(l-x- yf + (1 - yf - M 2 V I (2m\)\ 
(1 - 2y - m 2 /m' x ) 2 (3 - 2x - 2y + m 2 /ml) 2 ' 

(7) 



where x — Ey/m x , y — E q /m x , and x max = 1 + 
M£/(4m 2 ). The factor for the various channels reads 
C 7 = 3qja em , C z = 3q 2 f ta,n 2 (8 w )a em , and C g = 4a s , 
where qf — 2/3 for up- type and —1/3 for down- type 
quarks, a em — e 2 /(47r), and a s — g 2 /(47r). Since the 
gauge boson has to carry away one unit of angular mo- 
mentum in order to lift the helicity suppression, its spec- 
trum is rather hard and features a pronounced peak close 
to the dark matter mass [TjH [M] ■ When V is a photon, 
this process of internal bremsstrahlung can therefore lead 
to a characteristic feature in the gamma-ray spectrum 
that may be observable in current and future gamma-ray 
telescopes such as Fermi LAT or the planned Cherenkov 
Telescope Array (CTA) [94l [95] . In the present context, 
we concentrate on the antiproton flux, that is mainly 
produced by the annihilation into gluons, XX ~^ QQ9- 



^ r > 2m rLWi ¥ , (8) 

x 3 

where (cv}j is the thermally averaged cross-section mul- 
tiplied by the velocity in the annihilation channel j, p(r) 
is the distribution of dark matter particles in the Milky 
Way, and dNp/dT is the energy spectrum of antipro- 
tons produced in that channel per unit of kinetic energy. 
We assume for simplicity a spherically symmetric dark 
matter distribution, and calculate the antiproton flux as- 
suming a radial dependence given by either the isother- 
mal, Navarro-Frenk- White (NFW) [SSI [ST] or Einasto 
profiles [4Ql [98], with scale radii r s = 4.38, 24.42 and 
28.44 kpc [73j[74], respectively, and normalized to a local 
dark matter density p = p{r & ) = 0.4GeV/cm 3 [551 - 
[571 155] . The spectrum of antiprotons is obtained using 
the event generator CALCHEP [99] [TOO] interfaced with 
PYTHIA 8.1 [TOT] . 

The largest source of uncertainty arises from the prop- 
agation of antiprotons from the production point to the 
Earth. Following [102] , we use a stationary two-zone dif- 
fusion model with cylindrical boundary conditions. The 
number density of antiprotons per unit kinetic energy, 
fp(T, r, t), approximately satisfies the following transport 
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equation in this model: 



dfp 
dt 



V ■ (K(T,r)Vf p ) - V ■ (V c (r)f p ) 

-2hS(z)T ann f p + Q(T,r) . (9) 



The boundary conditions require the solution fp(T,r,t) 
to vanish at the boundary of the diffusion zone, which 
is approximated by a cylinder with half-height L = 1 — 
1 5 kpc and radius R — 20 kpc. The influence of the galac- 
tic magnetic field is modeled by an energy-dependent dif- 
fusion term, K(T) — Kof31Z s , where (3 — v/c is related 
to the velocity and 1Z — p(GeV)/\Z\ is the rigidity given 
by the ratio of momentum and electric charge. In addi- 
tion, the second term accounts for the drift induced by 
the Galactic wind, described by V c (r) = V c sign(z) k. 
Here k denotes the unit vector in the z-direction per- 
pendicular to the Galactic plane. The third term ac- 
counts for the annihilation of antiprotons when inter- 
acting with ordinary matter in the Galactic disk. The 
annihilation rate is taken from |73j . The flux arriving 
at the solar system is related to the number density by 
<&f{T) = £fp(T,r e ). Finally, the flux observed at the 
top of the Earth atmosphere is influenced by the effect of 
solar modulation, which we model within the force field 
approximation |103H105j . as detailed in [73]. Since the 
PAMELA data used in this analysis was taken near so- 
lar minimum activity, we choose tj>p = 500 MV for the 
solar modulation parameter for our numerical analysis. 
In order to take the uncertainties related to propagation 
into account, we use three sets of parameters, compati- 
ble with the cosmic boron-to-carbon flux ratio [1061 1107] , 
corresponding to minimum, medium and maximum an- 
tiproton flux with parameters as given e.g. in [73] , 



B. Experimental data &: background 

The most recent measurements of the antiproton 
flux by the PAMELA satellite experiment [37] in the 
energy range from 60 MeV to 180 GeV are in accor- 
dance with measurements performed by BESS [341 135] 
(0.1 - 4.2 GeV), IMAX [TUS] (0.25 - 3.2 GeV), WiZ- 
ard/CAPPJCE [109] (0.62 - 3.19 GeV), and AMS [HO] 
(0.2 - 4 GeV). In the future, data from the AMS-02 ex- 
periment are expected to yield further information on 
the cosmic antiproton flux (1111 1112j . The measured flux 
as well as the antiproton-to-proton ratio agree well with 
the expectations from secondary production of antipro- 
tons from spallation of cosmic-ray nuclei, mainly pro- 
tons and helium, on the interstellar medium |50l 1107] . 
This allows to set stringent upper limits on a possible 
primary contribution generated from dark matter anni- 
hilations. In our analysis we use the background flux 



calculated in Ref. [113] based on the two-zone diffusion 
model, taking into account p-p, He-p, p-He and He-He 
nuclear reactions. The main uncertainties arise from the 
diffusion parameters and the nuclear cross-sections, and 



are estimated to be in the range of 10 — 25% depend- 
ing on the energy. In contrast, the uncertainty stem- 
ming from the knowledge of the flux of cosmic nuclei and 
the composition of the interstellar medium are found to 
be subdominant. In order to obtain constraints on the 
dark matter annihilation cross-section, we compute the 
antiproton-to-proton ratio p/p = (Q p 6 + $£ ks )/$ p using 
the proton flux of [107j . In order to obtain a conserva- 
tive exclusion bound we adopt the minimal value for the 
antiproton background as discussed in (1 13] and a lower 
cut T > 1.5 GeV on the kinetic energy. Upper limits 
on the Yukawa coupling /, that controls the size of the 
annihilation cross-sections, are then obtained from the 
PAMELA p/p data [37] using a x 2 -test at 95% confidence 
level (CL). 



IV. DIRECT DETECTION 
A. Formalism & uncertainties 

The subject of direct detection of dark matter has been 
extensively covered in the literature [3J I114H116] , so here 
we limit ourselves to outline the features of interest for 
our study following closely |117[ I118j . Given a target 
nucleus N(A,Z), the differential rate of nuclear recoils 
induced by WIMPs is readily obtained by convoluting 
the local dark matter flux and the WIMP-nucleus scat- 
tering cross-section. We take the local WIMP velocity 
distribution to be a truncated Maxwellian with escape 
velocity v esc and circular velocity vq- The SI form factor 
is implemented as in [114j . while the SD one is modelled 
following [1 19] whenever possible or using the simplified 
prescription of [114] otherwise. The nuclear spin expec- 
tation values can be found in |120j for various nuclei (for 
nuclei absent in that reference we use |121j ). Finally, the 
remaining ingredients to specify in order to obtain the 
recoil rate are the effective SD and SI WIMP-nucleon 
couplings, a Pi n and / p , n . 

At the microscopic level, spin-dependent WIMP- 
nucleus scattering is induced by the exchange of a scalar 
in the s-channel or a Z boson in the i-channcl between 
the WIMP and a quark |115] , In the framework of the 
class of models introduced in Section [H] the only scalar 
at play is r\ and there is no WIMP coupling to the Z bo- 
son. Accordingly, the effective axial-vector WIMP-quark 
coupling reads 
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r- 



(10) 



8 m,2 — (m x + m q ) 
and the effective spin-dependent WIMP-proton coupling 



P ^ V2G F 

q—u,d,s 



(ii) 



and likewise for a n , where the parameters Aq( p ' n > set the 
spin content of the nucleon. 
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Spin-independent scattering, on the other hand, results 
from WIMP-quark and WIMP-gluon interactions via the 
exchange of scalars and/or Higgses at both tree and loop 
level 115 . Specialising to the particle physics realisa- 
tion in study - which features the scalar 77 only and no 
significant WIMP-Higgs coupling -, the effective spin- 
independent WIMP-proton coupling is 
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where /jT and /jfg, parametrise the quark and gluon 

content of the proton respectively, and q^(2), q( p \2) are 
the second moments of the parton distribution functions. 
A similar expression holds for /„. The effective WIMP 
couplings g q and b are given by 



9q 



(to 2 , - (m„ 



m q y 

™2 



B, 



-D 



2S 



B 1S cx f 



(13) 



in which B x are loop amplitudes that can be found in 
[1151 1122] . In Eq. (12 1, the first term corresponds to 



WIMP-quark scattering through the exchange of 77 in 
the s-channcl, the second term encloses the contribu- 
tion of WIMP-gluon scattering at one-loop level involv- 
ing quarks and 77, and the last term includes the twist-2 
operator part. The eventual contributions arising from 
couplings to the charm- and top-quarks are left out since 
they are absent of our scheme of couplings described in 
Section HU 

It is clear from the framework outlined above that the 
presence of a mass state rj almost degenerate with the 
dark matter particle strongly enhances both SD and SI 
direct detection rates. Indeed, the SD effective coupling 



d q and the SI term g q in Eqs. ( 10 ) and (Il3j) exhibit a 



resonance at — m x + m q . This was the motivation 
for the approach followed in [75] and is further pursued 
in the present work. In order to avoid modelling the 
width of the exchanged scalar rj - and thus introducing 
more freedom in our particle physics framework -, we 
restrict ourselves to mass splittings not too close to the 
resonance, namely (m^ — m x )/m q > 2 for each quark cou- 
pling considered. Also, we shall only derive constraints 
for — m x > 1 GeV ~ m p where the expansion lead- 
ing to the effective Lagrangian with couplings given by 
Eqs. ( 10 ) and ( 13 1 is valid. Even with such conservative 



restrictions, the amount of enhancement produced in di- 
rect detection by a compressed mass spectrum is very 
significant as shown in Section [Vj 

Before proceeding with setting constraints, it is in- 
structive to discuss the direct detection phenomenology 



of the dark matter models under consideration. First 
of all, we comment on the relative strength of each of 
the three SI terms in Eq. (12 1. Generically speaking, 



for the scheme of possible coupling configurations dis- 
cussed in Section [IT] the 6-term is always subdominant 
except for bottom-quark couplings. For the up-quark, 
down-quark and democratic light quarks couplings, the 
last term in Eq. ( 12 ) - i.e. the twist-2 contribution - 



plays the major role, while for strange-quark couplings 
it is the first term due to scalar exchange that domi- 
nates. Despite this sort of interplay, the scalar exchange 
and twist-2 contributions are roughly of the same order 
of magnitude, which makes direct detection constraints 
almost independent of the coupling scheme (for the case 
of light quarks). Secondly, the ratio of SD to SI WIMP- 
proton cross-sections a p /a p is easily much larger than 
unity, reaching ~ 10 2 — 10 s for the up-quark coupling with 
m rj jm x = 1.1 and using mean values for the nuclear pa- 
rameters (see below). Although the SD WIMP-nucleon 
cross-section is easily larger than its SI counterpart, it 
must be recalled that in direct detection SD sensitivities 
lag far behind SI searches. Actually, the strongest direct 
detection constraints on our DM model will be largely 
coming from SI couplings. Finally, let us point out that 
fiducial values are a n /a p = —0.5 and f n /fp = 0.6 in the 
case of up-quark couplings, and f n /fp = 1 in the case 
of bottom-quark couplings. These values depend on the 
adopted nuclear parameters, but are largely independent 
of and m x . 

Ultimately, our aim is to link the nuclear recoil event 
rate observed in a given experiment to the couplings and 
masses of the DM model. To pursue this path close atten- 
tion must be paid to two kinds of inputs: astrophysical 
quantities and nuclear parameters. On the astrophysical 
side, the key quantities are the local dark matter density 
Po and the local velocity distribution which is regulated 
by the circular and escape velocities, vq and v esc . We use 
the following fiducial values for these quantities |117] : 

po = 0.4 GcV/cm 3 , 
v Q = 230 ± 30 km/s , 
v esc = 544 km/s . (14) 

The largest astrophysical uncertainties under our frame- 
work regard po and vq. Since the former affects similarly 
direct detection and antiproton constraints, we opt to 
gauge its uncertainty out and stick to the typical value 
0.4 GeV/cm 3 [55H571 [55]. The uncertainty on v re- 
ported above is a reasonable assessment of the present 
(lack of) knowledge on this important galactic parame- 
ter [612H561ITTT] and shall be used throughout to represent 
the astrophysical uncertainty in direct detection. Other 
astrophysical parameters as v esc are also uncertain to a 
certain level [ST] , but such uncertainty affects recoil spec- 
tra much less and therefore shall be disregarded in the 
following. 

The input from nuclear physics amounts to the param- 
eters AqM in SD scattering and fft' n \ ffy n \ q& n )(2) 
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(P,n) 



and <p'™)(2) in the SI case. The coefficients 

and fj?G^ can furthermore be expressed in terms of the 
7T— nucleon sigma term E„.„ and the parameter <7 (see 
for instance )• The fiducial values of all these quanti- 
ties can be found in 1123] . Following [124j , the most 
important nuclear nuisances regard As' ? ' n ', E,™ and <To 
whose uncertainties we take from [69] : 

As (p) = As {n) = -0.09 ±0.03, 
= 64 ± 8 MeV , 
cr = 36 ±7 MeV . (15) 

Moreover, the parton distribution functions are particu- 
larly prone to nuisances and degrade the accuracy of SI 
searches. From Fig. 9 in |68j the relative uncertainty on 
the parton distribution function u(x) (d(x)) is of order 
10-25% (10-50%). Upon integration to get u{2) and 
<i(2), these uncertainties are partly smoothed out. Con- 
servatively, we consider that all the second moments in 



the last term of Eq. ( 12 1 are known to a 30% accuracy. In 



the remainder of the work, all these nuclear uncertainties 
are folded in with the astrophysical nuisances described 
in the previous paragraph to quantify the total uncer- 
tainty affecting direct detection. 



B. Experimental data 

In order to derive constraints on the parameter space of 
our DM model, we make use of the latest results from di- 
rect detection: XenonlOO [TH] for SI scattering, XenonlO 
[T7] for SD-neutron scattering, and SIMPLE [25] and 
COUPP [27] for SD-proton scattering. To the best of 
our knowledge these are the strongest published limits 
for each type of couplings. Recall that our focus here is 
on WIMP masses in excess of 40 GeV; data sets other 
than the ones mentioned above can of course be more 
constraining at lower masses. We briefly summarise be- 
low each experimental setup: 

• XenonlOO 19 . Using Xe as target material, the 
XenonlOO experiment has set the strongest direct 
detection limits to date by collecting 3 WIMP- 
like events in the energy window Er = 8.4 — 44.6 
keV for an effective exposure of 1471 kg. day. For 
such exposure the background estimate is 1.8±0.6 
events. Following the Feldman-Cousins procedure 
[125] with N bs = 3 observed events and Nbkg — 
1.8 mean expected background, we derive a 95% 
CL upper limit on the WIMP signal and obtain 
Nr < 6.45. Notice that, besides its superb sen- 
sitivity to SI couplings due to the high number of 
nucleons, xenon also presents good SD-neutron sen- 
sitivity through the abundant odd isotopes 129 Xe 
and 131 Xe. 

• XenonlO [T7]. Early results from the XenonlO 
experiment corresponding to an effective exposure 



of 136 kg.day show N obs = 10 events in the WIMP 
signal region within the energy range Er = 4.5 — 
27 keV. Since no background estimate is available, 
we stand on the conservative side and assume zero 
background, Nbkg — 0. Then, the Feldman-Cousins 
procedure yields an upper limit Nr < 17.82 at 95% 
CL. 

• SIMPLE [25J. The SIMPLE experiment consists 
of a threshold bubble device using C2CIF5 as tar- 
get material. We focus here on the results of phase 
II which comprises stages 1 and 2, with a com- 
bined after-cuts exposure of (13.47+6.71) kg.day. 
For this exposure, N a b s = 14 + 1 WIMP-like events 
were registered, while the neutron expected back- 
ground amounted to 0.896 counts/kg/day in stage 1 
and 0.253 counts/kg/day in stage 2 yielding Nbkg — 
12.07 + 1.70. Since the threshold recoil energy for 
bubble nucleation is E t hr = 8 keV in this experi- 
mental setup, we consider a range Er = 8 — 100 
keV for the derivation of our limits along with a 
bubble efficiency rf = 1 - cxp (— r {E R /E thr — 1)) 
where T — 4.2. Applying the Feldman-Cousins ap- 
proach with N obs = 15 and N bkg — 13.77 leads to 
an upper limit N R < 10.54 at 95% CL. This is a 
somewhat conservative limit because we do not ad- 
dress the problematics of double-scatterings, but it 
is sufficiently strong for our purposes. 

• COUPP [27]. Recently, the COUPP experiment 
- consisting of a CF3I bubble chamber - has re- 
leased results on WIMP searches for three distinct 
threshold energies, namely E t hr — 7.8, 11.0 and 
15.5 keV with total exposures 70.6, 88.5 and 394.0 
kg.day, respectively. The latter run turns out to be 
the most constraining for our purposes, and conse- 
quently we shall focus on such experimental setup 
where N obs — 8 WIMP-like events were observed. 
Following the analysis presented by the COUPP 
collaboration [27], a single-event efficiency of 79.1% 
is used to correct for the effective exposure and 
the bubble efficiency rf is set to unity for I re- 
coils above threshold and to 0.49 for C and F re- 
coils above threshold. Similarly to the case of SIM- 
PLE, an energy window spanning Er — 15.5 — 100 
keV is considered. Using the conservative hypoth- 
esis of zero background N b k g — 0, we employ the 
Feldman-Cousins procedure with N obs — 8 to de- 
rive N R < 15.29 at 95% CL. 



V. RESULTS & DISCUSSION 

We present in this section the results of our numerical 
analysis comparing the limits on the parameter space of 
our toy model which stem from the non-observation of 
an antiproton excess in the PAMELA cosmic antiproton- 
to-proton fraction or a signal in direct detection exper- 
iments. We show in Fig. [2] the limits on the coupling 
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FIG. 2: 95% CL constraints on / for the case of up-quark couplings and mass splittings m v /m x — 1.01,1.1,2. The yellow 
band encompassed by the dashed lines represents the antiproton constraints and corresponding uncertainties, while the direct 
detection limits are shown by the red band and solid lines. The thick dotted lines indicate the coupling / required to have a 
dark matter particle with the relic abundance measured by WMAP. Also shown as horizontal grey lines are the couplings of 
the supersymmetric bino to right-handed up-quarks, right-handed down-quarks and left-handed quarks. 



/ as a function of the dark matter mass coming from 
antiprotons (yellow band, dashed lines) and from direct 
detection (red band, solid lines), for various values of 
m ri /m x and assuming that the Standard Model fermion 
is a right-handed up-quark. The bands were constructed 
by marginalising over the uncertainties discussed in Sec- 
and 



tions 



III 



IV 



at each value of m x . In the case of 



matter mass m x > 45 GeV, 215 GeV and 830 GeV for 
m v /m x — 2, 1.1 and 1.01, respectively, as clear in Fig. [2] 
Notice that the quoted behaviour of the upper limits on 
/ with mr)/m x is rather strong if translated to cross- 
sections since a^ 1 oc f 4 . 

The limits on the parameter space of the model from 
antiprotons are generically weaker than those from di- 
rect dark matter searches. This can be explicitly seen in 
Fig.|3l whe re we plot the relative importance of each con- 
straint in the mass splitting vs dark matter mass plane. 
The red region delimits the parameter space for which di- 
rect detection limits take the lead over antiprotons even 
when uncertainties are taken into account - i.e. when 
the most conservative direct detection bound is stronger 
than the most aggressive antiproton constraint. The yel- 



direct detection, we chose the most constraining data set 
at each m x , which turned out to be always XenonlOO 
for the dark matter models studied here. Our numerical 
calculation shows that there exists very stringent limits 
on the coupling constant from direct detection experi- \3 
ments when the dark matter particle is close in mass to ^ 
the intermediate scalar particle. Whereas in the non- 
degenerate case, m v /m x = 2, the direct detection upper od 
limit on the coupling is / < 0.8 (10) for m x = 100 (1000) S 
GeV, in the degenerate case these limits get strength- 
ened to / < 0.08 (3) in the case m r] /m x = 1.1 and 
/ < 0.008 (0.4) in the case m v /m x = 1.01. As a compari- 
son, in the MSSM the Yukawa coupling of the bino to the 
up-quark and its corresponding squark is fsusY ~ 0.33, 
which would translate into a lower limit on the bino dark 




/// 



[GeV] 



FIG. 3: Most dominant constraints on the mass splitting vs 
dark matter mass plane for the case of up-quark couplings. 
The red (yellow) region signals the parameter space where 
direct detection is more (less) constraining than antiprotons. 
The grey region shows the region where direct detection and 
antiprotons are equally important within uncertainties. In 
the patch above and to the right of the dotted line thermal 
candidates present coupling constants in excess of 10. 



low regions refer to the opposite case, while the grey zone 
labelled "direct detection and antiprotons" marks the re- 
gion where antiprotons and direct detection are equally 
important within uncertainties. Only at large values of 
the dark matter mass the limits from antiprotons become 
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FIG. 4: The same as in Fig. [2] but for a mass splitting m^/m x = 1.1 and several coupling configurations. 



more important, due to the different scaling of the lim- 
its with the mass. However, with the current data, this 
occurs for values of the coupling constant where our per- 
turbative calculation is no longer valid - see dotted line 
in Fig. [3] Furthermore, as pointed out in Section |TV} we 
note that when the dark matter mass is close to the mass 
of r], the perturbative calculation of the direct detection 
rate breaks down, which happens at m x < 100 GeV for 
m n jm x = 1.01, hence in this regime only the antiproton 
limits are reliable. 

We also show in the plots of Fig. [2] as a dotted line, the 
coupling / required to give the dark matter thermal relic 
density as measured by WMAP. As discussed in Section 
HU this requirement sets a lower limit on the dark matter 
mass, which increases as m ri /m x approaches 1, namely 
m x > O[200 (1000) GeV] for m v /m x = 1.1 (1.01). Above 
the dotted line in Fig. [5] additional production mecha- 
nisms should be present, such as non-thermal production, 
in order to reproduce the observed dark matter relic den- 
sity. As can be seen from Fig. |2j in the degenerate case 
direct detection experiments start to probe the region 
of the parameter space favoured by thermal production 
and, in the non-degenerate case, even exclude the region 
m x < 100 GeV. 

To examine the sensitivity of these conclusions to the 
flavour of the quark to which dark matter couples, we 



have repeated the same analysis for couplings to the 
down-quark, the strange-quark, the bottom-quark and 
a democratic coupling to the up-, down- and strange- 
quarks; the results are shown in Fig. [4] for the case 
m ri /m x = 1.1. The conclusions for the three light quark 
flavours are very similar, the main difference being the 
width of the uncertainty band in the direct detection 
limits, which is larger in the case of the strange-quark. 
Furthermore, the antiproton and direct detection limits 
for the scenario where the dark matter particle couples 
democratically to the three light quark flavours are sim- 
ilar to those for the coupling to the individual flavours. 
Lastly, we have also calculated the limits on the coupling 
in a scenario where the dark matter particle only cou- 
ples to the bottom-quark. While the antiproton limits 
arc only slightly stronger than those for the couplings 
to the light quark flavours, the direct detection limits 
are significantly weaker. Still, the limits on this scenario 
from direct detection are stronger than from antiprotons. 
This is a particularly strong conclusion given the numer- 
ous works in the literature studying the complementarity 
between indirect and direct detection assuming WIMP 
annihilations into pairs of bottom-quarks. 

The previous limits on the coupling / can also be 
translated into upper limits on a^ 1 and (crv). The for- 
mer translation suffers, of course, from the nuclear un- 
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FIG. 5: Antiproton and direct detection 95% CL constraints on the total annihilation cross-section (left) and on the spin- 
independent WIMP-proton scattering cross-section (right) for the case of up-quark (top) and bottom-quark (bottom) couplings 
and a mass splitting m v /m x = 1.1. The colour and line code is the same as in Fig. [5] 



certainties discussed in Section IIVI The limits from di- 
rect detection and antiprotons are presented in Fig. [5] 
as a red and a yellow band respectively, for the case 
m r] /m x = 1.1 and couplings to up- and bottom-quarks. 
As mentioned above, for low dark matter masses direct 
detection experiments provide the most stringent con- 
straints on the model, while for large masses antiproton 
searches become dominant. However, we stress that the 
region at large masses where the antiproton limits are 
stronger than the direct detection limits corresponds to 



couplings larger than 10, hence the perturbative calcula- 
tion is no longer reliable. Moreover, it should be borne in 
mind that the bands shown in Fig. [5] are sensitive to as- 
trophysical quantities: the antiproton limits can become 
stronger in the presence of "boost factors" and the direct 
detection limits can become stronger or weaker depend- 
ing on the concrete value of the dark matter density at 
the solar system. We finally notice that the gamma-ray 
and antiproton limits on annihilations into quarks are 
of similar strength [95 , and therefore any interpretation 
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of the gamma-ray feature around 130 GeV as internal 
bremsstrahlung \X ~^ 997 1S m tension with direct de- 
tection constraints. 




m x [GeV] 

FIG. 6: 95% CL exclusion regions on the mass splitting vs 
dark matter mass plane for the case of up-quark couplings 
and fixing / = 1. The solid red, thin dashed yellow and thick 
dashed green contours correspond to the regions of the pa- 
rameter space excluded by direct detection, antiprotons and 
collider searches, respectively. The upper solid contour shows 
the reach of sensitivity of XenonlT. 

We now turn to the complementarity of antiproton and 
direct detection constraints in light of collider searches. 
For concreteness, we shall focus on the case of a WIMP 
coupling to up-quarks only. Fig. [6] shows the regions in 
the mass splitting vs dark matter mass plane excluded by 
direct detection (red, solid) and antiprotons (yellow, thin 
dashed) assuming a coupling / = 1 and taking the most 
conservative limits - i.e. the upper lines in the bands of 
Fig. [2j At m x = 300 (1000) GeV direct detection ex- 
cludes mass splittings m v /m x — 1 < 19 (2)% for / = 1 
and 6 (< 1)% for / = fsusv = 0.33. Notice that the 
direct detection contour towards the bottom left corner 
of Fig. [6] signals the limit of validity of our calculations 
(m v — m x > 1 GeV); the region below this line, al- 
though not strictly excluded by direct detection in our 
analysis, should be regarded as very contrived. Over- 
plotted to the antiprotons and direct detection limits in 
Fig. [6] is the region excluded by ATLAS (green, thick 
dashed) in the search for jets plus missing transverse mo- 
mentum |126| . Additionally, a recent re-analysis |127| of 
monophoton plus missing transverse energy data derives 
a lower limit m n > 150 GeV for up-type quarks provided 
m„ — m x < 10 GeV. These constraints, however, only 
probe a small patch in the bottom left corner of Fig. [6] 



which is largely excluded by direct detection if / = 1. 
Fig. [6] clearly shows that collider searches are effectively 
orthogonal to antiproton and direct detection limits. In 
fact, although the allowed region is large for high masses, 
in the low-mass regime only a small region still survives. 
This underlines the importance of pursuing complemen- 
tary searches to close in on mass-degenerate dark matter 
models. 

It is therefore interesting to investigate the prospects 
for antiproton, direct and collider searches. We start 
with antiprotons. The AMS-02 detector [T2"g] - installed 
on the International Space Station in 2011 - is presently 
collecting data and will likely extend the measurement 
of the antiproton flux to higher energies than PAMELA. 
However, it will probably take some time to reach the 
accuracy of PAMELA antiproton-to-proton data, which 
is the data we employ in the present analysis. Therefore 
we assume that antiproton constraints on dark matter 
models will not improve significantly in the short term. 
Direct detection, instead, presents promising prospects 
with several ton-scale experiments |129[ I130j already on 
the way. We focus on a relatively short time scale and 
take XenonlT as an example: it is expected [131] that 
XenonlT will improve a factor ~ 200 in SI cross-section 
sensitivity with respect to present XenonlOO results. As 
shown in Fig. |6j such prospects will effectively shift the 
direct detection constraints (which are dominated by 
XenonlOO) upwards and, assuming a coupling / = 1, 
these will cover the whole reach of current LHC searches. 
It is clear that a big chunk of the parameter space of 
mass-degenerate models will be excluded over the next 
few years if no signal is detected. In that case we expect 
stringent lower limits on the mass splitting of order 114% 
(10%) at m x = 300 (1000) GeV for / = 1. 



VI. CONCLUSION 

We have shown in this work that the latest data on 
cosmic antiprotons and direct detection place useful con- 
straints on the phenomenology of mass-degenerate dark 
matter scenarios. In particular, we have considered a 
minimal framework featuring a Majorana fermion as dark 
matter that couples to light quarks via a scalar close in 
mass, encompassing e.g. a simplified model with bino- 
like neutralino and squark as lightest and next-to-lightest 
supersymmetric particles. This setup allows for a di- 
rect comparison of scattering rates on nuclei with dark 
matter annihilation rates in our Galaxy, the dominant 
channel being \X ~~ ^ QQ9 f° r quarks of the first and sec- 
ond generation. The derived constraints on coupling to 
quarks suffer from sizeable astrophysical and nuclear un- 
certainties, but it is nevertheless clear that antiprotons 
lag significantly behind direct detection, a fact that can 
be attributed mainly to the extreme sensitivity of un- 
derground searches to mass degeneracy. Fine degenera- 
cies are conclusively discarded by current direct detection 
data. This is precisely the range that escapes detection 
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at collider searches. Accordingly, we find that the in- 
terplay between antiprotons, direct and collider searches 
will be of crucial importance in closing in on simple mass- 
degenerate dark matter models over the coming years. 
Further work is needed to study the implications of this 
complementarity in the framework of more complicated 
particle physics realisations. 
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